Atomic clock transitions in silicon-based spin qubits 
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A major challenge in using spins in the solid state for quantum technologies is protecting them 
from sources of decoherence. This can be addressed, to varying degrees, by improving material 
purity or isotopic composition [31^ for example, or active error correction methods such as dynamic 
decoupling [3114], or even combinations of the two [5l[6]. However, a powerful method applied to 
trapped ions in the context of frequency standards and atomic clocks [7| E] , is the use of particular 
spin transitions which are inherently robust to external perturbations. Here we show that such 'clock 
transitions' (CTs) can be observed for electron spins in the solid state, in particular using bismuth 
donors in silicon [5] [TO] . This leads to dramatic enhancements in the electron spin coherence time, 
exceeding seconds. We find that electron spin qubits based on CTs become less sensitive to the local 
magnetic environment, including the presence of ^^Si nuclear spins as found in natural silicon. We 
expect the use of such CTs will be of additional importance for donor spins in future devices 
mitigating the effects of magnetic or electric field noise arising from nearby interfaces. 



Out of the various candidates for solid state qubits, 
spins have been of particular interest due to their rel- 
ative robustness to decoherence compared to other de- 
grees of freedom such as charge. So far, the most coher- 
ent solid state systems investigated have been the spins 
of viTcU-isolated donors in bulk 28-silicon, with coherence 
times (T2) of up to seconds (extrapolated) for the electron 
spin [T] and minutes for the nuclear spin [5], comparable 
to those of ion trap qubits OII^. However, in practical 
devices, spin coherence times are likely to be limited by 
factors such as coupling to nearby qubits and magnetic 
or electric field noise from the environment. For exam- 
ple, cross-talk with other donors 100 nm away limits the 
electron spin to a few milliseconds [1 , while a nearby 
interface can limit the donor electron spin to 0.3 ms 
at 5.2 K p/P. Finally, without isotopic enrichment, the 
5% natural abundance of ^^Si limits the electron spin 
T2c to less than 1 ms [IHU]. 

An approach to creating more robust qubits is to tune 
free parameters of the system Hamiltonian to obtain in- 
sensitivity to specific sources of decoherence. This has 
been extensively used in ion trap qubits to protect against 
magnetic field fiuctuations [121 113j , building on work on 
atomic clocks where hyperfine states, used as frequency 
standards, must remain stable against such variations. 
These so-called "clock transitions" (CTs) have a transi- 
tion frequency (/) which is insensitive to magnetic field 
(B) variations, at least to first-order (in other words 
df/dB = 0). More recently, superconducting circuit 
qubits have also taken advantage of a tuned Hamiltonian 
to remain immune to charge, fiux or current noise |15U16j . 



Nuclear spin CTs in rare-earth dopants (nuclear spins 
/ > 5 /2) have been studied in the context of optical quan- 
tum memories [T71 [TS] leading to a 600-fold improvement 
of the coherence times to 150 ms, limited by second-order 
effects, while recent experiments on phosphorus donor 
nuclear spins also exploited a CT [5]. For electron spins 
in the solid-state, CTs remain relatively unused due in 
part to the requirement of a spin Hamiltonian of sufficient 
complexity. One of the richest single-defect spin systems 
is the bismuth donor in silicon (Si:Bi), which possesses an 
electron spin 5 = 1/2 coupled to a nuclear spin / — 9/2. 
The electron spin decoherence rates for Si:Bi have been 
found to follow df/dB in both natural silicon and 
isotopically enriched ^^Si [15]. These results, combined 
with the identification of a number of CTs in the spin 
Hamiltonian of Si:Bi ^20,i21J, motivate the study of spin 
coherence times around CTs in Si:Bi, where df/dB — > 0. 
In this Letter, we investigate one such CT in Si:Bi, at 
7.0317 GHz, using both natural silicon and ^^Si. 

When describing the states of coupled electron and nu- 
clear spins, two basis conventions are typically used: in 
the high magnetic field limit, the electron and nuclear 
spin projections ms and mj are good quantum numbers, 
while in the zero-field limit, the total spin F (= I zL S) 
and its projection mp {— ms + mj) are used. CTs are 
often found in an intermediate regime |22| . nevertheless 
it is possible to categorize them as nuclear magnetic reso- 
nance (NMR)- or electron spin resonance (ESR)-type, on 
the basis of whether the transition couples primarily to 
Sx or Ix , where these are the electron and nuclear spin op- 
erators perpendicular to the applied magnetic field. The 
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FIG. 1. Electron spin resonance (ESR)-type clock transitions (CTs) of Si:Bi. A, The eigenstate energies (top) of 
Si:Bi as function of magnetic field, tiie ESR- (black) and NMR-type (grey) transition frequencies between these states (middle), 
and the first-order magnetic field dependence (df/dB) of these transition frequencies (bottom). ESR- type CTs (blue lines and 
open circles) are found at 27, 80, 133 and 188 mT, and appear in the spectrum as doublets AFAmp — ±1 separated by up to 
3 MHz. NMR-type CTs are found above 300 mT (red lines and open circles). B, Electron spin echo-detected magnetic field 
sweeps around the 80 mT CT measured at microwave frequencies > 7.0315 GHz. The transition probablities for AFAmp = -1-1 
(dark blue) and —1 (light blue) transitions are equal near the CT. 



ESR-type CTs which we investigate in this manuscript 
involve states which are close to pure in the m^?) ba- 
sis and hence for convenience we label them according to 
the dominant \F,mp) component (full details are given 
in the Supplementary Material and in Ref |20j). 

For bismuth donors in silicon, NMR-type CTs can 
be found at high field (> 350 mT) with frequencies 
around 1 GHz as shown in red in Figure [ipL. At low 
field (< 200 mT), four ESR-type CTs are present with 
frequencies in the range 5.2 to 7.3 GHz as shown in blue 
in the same figure. We will focus here on the ESR-type 
CTs, which possess only slightly reduced spin manipu- 
lation time compared to free electron spins as well as a 
large energy splitting even at low magnetic field (which 
has interesting applications for use in hybrid supercon- 
ducting circuits 

In the silicon samples we study here, Bi donors were 
introduced during crystal growth using the method de- 
veloped in Ref [25 , with concentrations ranging from 
3.6 X 10^4 cm-3 to 4.4 x 10^^ cm^^. Pulsed-ESR experi- 
ments were performed using a spectrometer based around 
a modified Bruker Elexsys E580 system with a ~7 GHz 
loop-gap cavity (for the CT) and 9.75 GHz dielectric res- 
onator. 

Figure [lj3 shows ESR spectra measured using mi- 
crowave frequencies between 7.031 and 7.034 GHz, by 
plotting electron spin echo intensity as a function of mag- 
netic field. The spectra show two transitions correspond- 



ing to [{AF,AmF} = {±1,±1}] and [{Ai^,AmF} = 
{±l,=pl}]; for brevity, these transitions can be distin- 
guished by the value of the product AFAmp = ±1. 
Together, they offer a controllable two-qubit subsystem 
with low sensitivity to magnetic field fluctuations (see 
inset of Figure [Ij3). 

We model the ESR spectra using an isotropic spin 
Hamiltonian common for group V donors in silicon: 

Ho = Boi-/eS,®t--fnl®I,)+AS.I (1) 

where the two first terms correspond to the electronic 
(S) and nuclear (/) spin Zeeman interactions with an 
external field Bq and the last term corresponds to the 
hyperfine coupling A. A common way to estimate Hamil- 
tonian parameters such as the electron and nuclear gyro- 
magnetic ratios (7e and 7„) and the hyperfine constant is 
by measuring the magnetic field dependences of the spin 
transition frequencies. We use the opportunity provided 
by the CT (with df/dB — >■ 0) to extract a measure of the 
hyperfine constant A = 1.47517(6) GHz with high preci- 
sion, because uncertainties in the magnetic field become 
irrelevant. In our simulations, we additionally use the 
previously reported value of 7^ = 27.997(1) GHz/T ^ 
and the generic value of 7„ = 7 MHz/T for ^osBi [9]. 

Figure [lj3 shows that the ESR linewidth in the mag- 
netic field domain increases around the CT: the deriva- 
tive df/dB tends to zero hence its inverse, dB/df, di- 
verges until it becomes limited by the non-linear terms 
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FIG. 2. Decoherence mechanisms of Bi donors in sil- 
icon and their dependence on df/dB. A, In the central 
spin representation, a Bi donor is coupled to neighbouring Bi 
donors as well as ^®Si spins. At the ESR CT, all spectral dif- 
fusion (SD) contributions to decoherence are essentially elim- 
inated, leaving only the direct flip-flop term (dFF) between 
the central spin and a neighbouring, resonant Bi spin. B, T2c 
measurements at 4.8 K show a strong dependence on df/dB, 
as shown for 3 different donor concentrations in ■^®Si:Bi. Mea- 
surements close to df/dB — 7e were taken using the ten X- 
band ESR transitions, while the remaining points were taken 
close to the CT. For each concentration, the dependence on 
df/dB is modeled using contributions from ID, FF and IFF, 
as shown separately in dashed lines for the lowest concen- 
tration. Inset shows the limit of I/T2C when approaching to 
the exact CT as a function of donor concentration, showing a 
nearly linear dependence, as expected for dFF. 



in /(i?o)- These spectra are all well fit assuming a con- 
stant linewidth in the frequency domain of 270 kHz. This 
hnewidth can be attributed to a distribution in the hy- 
perfine constant of around 60 kHz, using A/ = ^AA 
at the CT. Fourier- Transform ESR performed at a range 
of frequencies confirmed that the ESR hnewidth in fre- 
quency domain is indeed magnetic field independent (see 
Supplementary Material). 

We now examine the decoherence mechanisms which 
affect the electron spin of donors in silicon. At sufficiently 



low temperature (< 5 K), spin-lattice relaxation Tic can 
be mostly neglected, and dipolar interactions (~ 2SzSz — 
(SxSx + SySy)) with neighbouring spins are the primary 
source of decoherence. In a central spin representation, 
as shown in Figure [2]\, the surrounding spins can be 
divided into three categories: i) resonant spins affected 
by microwave excitation; ii) off-resonant spins of the same 
species, i.e. Bi spins in mp levels not addressed by the 
microwaves; and iii) other spin species such as ^^Si. Away 
from CTs, the limiting factor for electron spin coherence 
times is spectral diffusion (SD) from the SzSz term of the 
dipolar interaction. This term can be assimilated into 
effective fluctuations in the magnetic field environment 
of the central spin. SD is independent of any frequency 
detuning between spins and thus is valid between the 
central spin and any others. 

In the static case, dipolar couplings to (ii) and (iii) 
can be refocused with a microwave 7r-pulse such as in 
the Hahn echo sequence. However, this does not cor- 
rect for the dipolar coupling between resonant spins (i) 
as both spins are simultaneously flipped by the 7r-pulse. 
This is called "instantaneous diffusion" (ID) and lim- 
its T20 to ~ 10 — 100 ms for typical donor concentra- 
tions (> 10^'' cm~'^) [HITn] [37]. Furthermore, dynamic 
changes from spin flips in the environment cannot be refo- 
cused. At high temperature, such flips arise from phonon 
scattering but at low temperature, this is due to flip-flops 
(FF) from the SxSx + SySy term of the dipolar interac- 
tion. FF are energy conserving and as such are only rel- 
evant between spins that have similar transition frequen- 
cies. In natural silicon, the dominant decoherence mech- 
anism is SD from ^^Si FF, while in isotopically enriched 
^^Si, it arises from FF between resonant Bi spin pairs. In 
the latter case, we distinguish between FF which involve 
the central spin (direct FF, dFF), and those which do 
not (indirect FF, IFF). 

We begin by discussing results on samples of isotopi- 
cally enriched ^^Si (100 ppm ^^Si). At the CT the transi- 
tion frequency is insensitive to magnetic field fluctuations 
in first order, so we expect SD to have little effect, leav- 
ing only the dipolar coupling between resonant spin pairs. 
With reference to Figure [2]A, this implies then that all 
terms apart from dFF vanish. In Figure [2^3, measure- 
ments of electron spin coherence times (Tjjo) are shown 
for three different concentrations over a wide range of 
df/dB. The data includes values measured at X-band as 
weU as those near the CT {AFAmp = +1) at 79.8 mT, 
7.0317 GHz. Measurements at the CT shown here were 
taken at 4.8 K where Tic = 9 s, however no increase in 
was seen at lower temperature. 

For each sample, enhancements of about two orders of 
magnitude are seen at the CT, compared to the case for a 
free electron g-factor, such as that of phosphorus donors. 
As shown in Figure [2j3, the dependence of the measured 
on df/dB arises from two factors: the effect on ID, 
and on IFF. ID has a known quadratic dependence on the 



4 




2 4 
Time, 2r(s 



0.1 0.2 0.3 
Time, 2r(s) 



FIG. 3. Hahn echo decay at the CT. A, ^^SiiBi at 4.3 K 
with a Bi concentration of 3.6 x lO" cm"^. B, "^*Si:Bi at 
4.8 K with a Bi concentration of 10^^ cm~^. The decay in 
natural Si is a stretched exponential, and therefore T2c is de- 
fined as the time when the amphtude reaches 1/e. Magnitude 
detection was used to ehminate instrumental noise, likely due 
to phase noise in the microwave source. 



gyromagnetic ratio of the central spin [THl , and be- 
comes a negligible effect for df /dB < O.l^e- Indirect FF 
dephase the central spin through the SzSz term, giving 
a linear dependence of on df /dB. Direct FF, on the 
other hand, are not eliminated at the CT, and provide 
an upper bound on T20 for a given donor spin concentra- 
tion, as plotted in the inset of Figure [2j3. For the lowest 
concentration sample, electron spin coherence times of 
up to 2.7 s were measured from simple two-pulse Hahn 
echo decays, as shown in Figure |3j^. 

We now turn to measurements on Bi-doped natural sil- 
icon ("^'Si:Bi), which has 5% ^^Si. Away from the CT 
the effect of the ^^Si (J = 1/2) is both to broaden the 
ESR linewidth to about 0.4 mT (equivalent to 12 MHz 
in the frequency domain for a free electron) due to unre- 
solved ^^Si hyperfine, as well as to limit the to about 
0.8 ms due to SD [5]. At the CT we find that the ESR 
linewidth reduces to 500 kHz (see Supplementary Ma- 
terial), within a factor of two of the value for enriched 
^^Si material, while increases by over two orders of 
magnitude to about 90 ms (Figure [sJB). The effect of 
the suppression of SD around the CT has been simu- 
lated for '^'^*Si:Bi using cluster expansion methods 
though further refinements are required in the simula- 
tion before a quantitative comparison can be made. The 
stretched-exponential decay implies that is still lim- 
ited at 93 ms by SD from ^^Si due to the second order 
term [d^f/dB^ ^ 0). For modest ^®Si enrichment (e.g. 
P^Si] « 1000 ppm), T2c should already exceed seconds, 
and indeed there may be an optimal ^*Si purity above 
which T2c at the CT drops, due to the role of ^^Si or ^°Si 
in detuning otherwise-identical spins |30j . 

We have shown how CTs in Si:Bi can be used to pro- 
duce magnetic field-insensitive spin qubits with directly 
measured coherence times of several seconds. Such qubits 
would be insensitive to magnetic field noise arising, for 
example, from fluctuating dangling-bond spins at the 



Si/Si02 interface. Conversely, if electric field noise is 
dominant, this can couple to donor spins via the hyper- 
fine interaction and cause decoherence. Again, CTs can 
be designed to be immune from electric charge noise by 
selecting points where df/dA — ^ (see Supplementary 
Material). Through the use of CTs, it is likely that the 
seconds-long electron spin coherence times measured in 
the bulk can be harnessed for spins in practical quantum 
devices. 
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ESR- AND NMR-TYPE MAGNETIC FIELD "CLOCK" TRANSITIONS (CT) 




Magnetic field (mT) 

FIG. 1. Description of ESR-type CTs. A, The eigenstate energies of Si:Bi as function of 
magnetic field. Tlie color scale shows the logarithmic distance to pure Bell states (Landau- 
Zener (L-Z) anticrossings) in the \ms,mj) basis, defined as log(|^ — it/4\) for an eigenstate 
|<I>) = cos{9) |=p^,m/ lb ^) lb sin(0) |ib^,m/ T 5). The Bell state at mT is barely visible here 
due to degeneracy. B, ESR-type CTs with AFAmp = +1 in dark blue and —1 in light blue. The 
four involved eigenstates, which are two pairs of hyperfine coupled states in the |m5,m/) basis, 
form a subspace of the Hilbert space. 

Below we discuss the general requirements for ESR and NMR-type CTs for systems with 
electron spin S = 1/2 and nuclear spin I and assuming an isotropic hyperfine coupling, 
with a particular focus on Group V donors in silicon. This letter is the first measurement of 
ESR-type CTs to our knowledge, though they have been theoretically described by Moham- 
mady et.al. [1, 2] for donors in Si, in particular Bi. On the other hand, NMR-type CTs have 
been used in various systems in the past [3, 4], including in phosphorus donors in silicon [5], 
to reduce sensitivity to magnetic field inhomogeneities (i.e. increase frequency resolution) 
or to increase nuclear coherence times. 

In the basis of the electron and nuclear spin ImsyTrij), the isotropic hyperfine interaction 
{AI ■ S) couples pairs of states within the Hilbert space such that [Ams = ±1, Am/ = =i=l]. 
In the strongly coupled electron-nuclear spin basis \F,mF) {F = I±S,mF = ms+mj), these 
pairs of states share the same mp value. When the static magnetic field is increased, the 
Zeeman energy rises to the same order of magnitude as the hyperfine interaction, resulting 

1 



in avoided Landau-Zener crossings between states with < as shown in Figure lA. 
ESR-type CTs are located between pairs of these avoided crossings. 

At the intermediate fields relevant to CTs and Landau-Zener anticrossing, the eigenstates 
can be expressed either as being close to Bell states in the |m5, m/) basis, or still quite pure 
in the |-F, mp) basis. For example, the CT at 7.0317 GHz, explored in the main text, 
connects the following pairs of states: 
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in the |ms,mi) basis 
0.99 |4, -2) + 0.15 |5, -2) ^ -0.15 |4, -1) + 0.99 |5, -1) in the |F, mp) basis 



Hence, for convenience, we refer to these states by the dominant term in the \F,mF) basis 
(i.e. |4, -2) and |5, -1) in the example above). The AFAmp = +1 and AFAmp = -1 CTs 
are each transitions between one state of the first Landau-Zener crossing to a second state 
of the second crossing, forming a 4-dimensional subspace of the Hilbert space (Figure IB). 

For non-integer [integer] values of /, there are / -1-1/2 [/] Landau-Zener anticrossings and 
consequently 2(J — 1/2) [2/] ESR-type CTs. The minimum complexity required is thus a 
nuclear spin / > 1 to have at least two pairs of hyperfine coupled states, which is not the 
case for phosphorus donors (^^P has 1 — 1/2). 

Arsenic (^^As, I = 3/2) and antimony (^^igb, I = 5/2 and ^^^Sb, I = 7/2) have sufficient 
nuclear spin to permit ESR-type CTs, however the hyperfine coupling is relatively weak 
(^4 ~ 198, 186 and 101 MHz, respectively), so the avoided crossings are found at low 
magnetic field and transition frequencies (see Table I). Bismuth ^"^Bi with 1 = 9/2 and 
A = 1.475 GHz is thus the optimal Group V donor in silicon from the point of view of CTs, 
possessing four of them at GHz frequencies. 





''8 As (/ = 3/2) 


i2iSb (7 = 5/2) 


i23Sb (7 = 7/2) 


209Bi (7 = 9/2) 


AF = +l,mF = 
Magnetic field (mT) 
Frequency (GHz) 


-1 •(-)■ 
3.8 
0.384 


-1^0 
3.4 
0.552 


-2 ^ -1 
10.4 
0.482 


-1^0 
1.8 
0.403 


-2 ^ -1 
5.5 
0.376 


-3 ^ -2 
9.3 
0.314 


-1^0 
26.6 
7.338 


-2 ^ -1 
79.8 
7.032 


-3 ^ -2 
133.3 
6.372 


-4 ^ -3 
187.8 
5.214 



TABLE I. Summary of ESR-type magnetic-field CTs in donors in silicon. CTs exist in 
pairs {AFAmp = ±1) separated by less than 0.15 mT in magnetic field and 3 MHz in frequency. 
Phosphorus does not possess any CT due to its small nuclear spin (7 = 1/2). 



As df /dB — > 0, the next figure of merit is the electron transition probability amplitude. 
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which is always 50% of the high field limit for ESR-types. This means that manipulation 
times are only slightly reduced while electron coherence times ^-re drastically increased. 
Conversely, as flip-flops between two donor electron spins follow the square of this probability 
amplitude, the flip-flopping rate remains strong, limiting as observed and explained in 
the main text. 

NMR-type CTs occur at strong magnetic held (in Si:Bi, 0.3 T < Bq < b T, from mi — 
—9/2 at low field to mj = 9/2 at high field) and as such are transitions between quasi- 
pure states in the (m^, mi) basis. NMR-type CTs possess a change in nuclear spin state 
of Ami — 1 and can be manipulated in a conventional electron nuclear double resonance 
(ENDOR) or NMR experiment. Amongst the NMR-type CTs, those at higher magnetic 
fields have a smaller electron spin component. This leads to a reduced coupling to the 
environment, but also increased spin manipulation times, converging to that of a regular 
NMR transition (typically ~ 10 //s). 



ESR LINEWIDTHS 



Measurements of spin linewidths provide important details regarding the spin environ- 
ment, yielding information on crystalline defects and strains amongst other properties. Line 
broadening arises from a variation in either the hyperfine interaction (AA) or the magnetic 
field {AB) (or g-tensor) across the sample. They can either be measured in a magnetic 
field-swept spectrum where (to first order): 

or in Fourier- Transform (FT) ESR where the FT of the free induction decay after a 7r/2 
rotation gives the spectrum in the frequency domain. In this case, (to first order): 

At the CT where dB/df — > oo, the linewidth in the magnetic field domain broadens 
strongly until it becomes limited by the second order dependence on magnetic field. Both 
AB and AA can be identified by varying the transition frequency and fitting numerically 
knowing ^ To confirm that this change in hnewidth is strictly related to the term df/dB, 
we can make use of the FT ESR technique which simphfies at the CT to Af = ^AA, where 




Frequency (MHz) Frequency (MHz) 

FIG. 2. FT-ESR around the CT. Each spectrum is the FT of the free induction decay taken 
using a microwave frequency shghtly below resonance (values given in legend). A, For the case of 
^*Si:Bi, we observe two peaks in the ESR spectrum around the CT, corresponding to the transitions 
l^F^mp = ±1- Spectra are shown as measured at two settings of magnetic field/microwave 
frequency. In the magnetic field domain, the ESR linewidths in these two cases are 1.6 mT close to 
the CT and 0.07 mT farther away (see Figure 1 of the main manuscript), however in the frequency 
domain as shown above, the ESR linewidths are constant. B, In "^*Si:Bi, these two primary ESR 
transitions are further split into sub-peaks, corresponding to a mass-effect from nearest neighbour 
Si atoms. Each shift of one neutron mass (A/i) yields a shift of —1.7 MHz in transition frequency 
(or 0.024% change in the hyperfine coupling A). Dashed lines show simulated peaks whose intensity 
is calculated from a trinomial distribution of and ^®Si isotopes in '^'^^Si (with respective 

concentration 3.1%, 4.7% and 92%). As the FT-ESR is derived from the free induction decay, the 
intensities are normalized by the FT of the inhomogeneous decay T2e* (Lorentzian) and the cavity 
bandwidth. 

df /dA is quasi-constant around the CT. In Figure 2A, the linewidth is indeed constant about 
270 kHz. 

In natural silicon, the elimination of the Ai? term dramatically reduces the ESR 
linewidth. Away form the CT (e.g. at X-band), AB is normally around 4 G due to unre- 
solved coupling to ^^Si nuclear spins. In the frequency domain, this would be equivalent to 
nearly 12 MHz, hiding multiple spectral features. First, the two transitions AFAmp = ±1 
would not be resolvable. Second, as show in Figure 2B, we observe several other peaks (ab- 
sent in isotopically pure ^^Si samples) which arise from variations in the hyperfine coupling 
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due to the total mass of nearest-neighbour sihcon atoms (^^Si, ^^Si and ^"Si). This effect is 
described in full detail, including ENDOR experiments, in a forthcoming work [6] . 

ELECTRIC FIELD CLOCK TRANSITIONS 




200 400 600 800 1000 
Magnetic field (mT) 



FIG. 3. Clock transitions in Si:Bi where df /dA = which should be robust to electric 
field noise. Both ESR- and NMR-type CTs can be observed, as for magnetic-field CTs. One 
further ESR-CT is found at higher magnetic fields (2.6 T, not shown). 

In the experiments reported in the main text, CTs were used to reduce the sensitivity 
of electron spin to magnetic field variations, as quantified by df/dB. While magnetic field 
noise is indeed the main decoherence mechanism in bulk materials, this may not be the case 
in nanoscale devices where the electric field at interfaces could couple strongly with both 
the donor electron and nuclear spins through the hyperfine interaction (and also, to lesser 
extent, through a modulation in the electron spin g-factor). The sensitivity of a spin to 
this effect can be quantified by the gradient of the frequency with respect to the hyperfine 
constant df/dA, combined with values for the DC Stark effect for donors in silicon (which 
for Group V donors is in the order of 10~^ yum^/V^, as a fractional change in the hyperfine 
coupling [7, 8]). Those CTs which will be most robust to electric field noise {df/dA — )■ 0) 
are identified in Si:Bi in Figure 3 and in Table II for all Group V donors in silicon. 

In practice, both magnetic and electric field fluctuations will participate to the donor spin 
decoherence. There will thus be an optimal CT, at a specific magnetic field and frequency, 
where the coherence time would be maximum. For example, the magnetic field CT near 
188 mT in Bi has the lowest value of df/dA out of the four possible CTs. In other scenarios, 
it might be advantageous to minimise the inhomogeneous broadening as much as possible 
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75 As (/ = 3/2) 


i2iSb (/ = 5/2) 


123Sb (7 = 


7/2) 


209Bi (7 = 9/2) 


Ams — +1, mi = 


-1/2 


-1/2 


-3/2 


-1/2 


-3/2 


-5/2 


-1/2 


-3/2 


-5/2 


-7/2 


Magnetic field (mT) 


53 


117 


39 


114 


38 


23 


2607 


868 


519 


369 


Frequency (GHz) 


1.43 


3.21 


0.92 


3.17 


0.98 


0.49 


72.64 


23.18 


12.57 


7.30 



TABLE II. Summary of ESR-type electric-field CTs in donors in silicon. At the given 

magnetic fields, the electron and nuclear spins are weakly coupled and the eigenstates must thus 
be expressed in the \ms,mi) basis. The [Ams = ±l,Am/ = ^2] {AFAmp = —1) transitions 
are nearly completely forbidden here; they would have been found at the same magnetic field as 
the [Ams = ±1, Am/ = 0] {AFAmp = +1) transitions, but separated by less than 40 MHz in 
frequency. 

(e.g. for coupling a spin ensemble to a microwave resonator), and this would also require 
different optimal operating points within the Hilbert space of the bismuth electron and 
nuclear spins. 
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